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approach. The quark correlator is first investigated in the framework of QCD sum 
rules and the correlation length of the quark field is calculated. Comparing the 
phenomenological part of the sum rule with previous measurements of the quark 
correlator on the lattice, we are able to obtain an independent result for the corre- 
lation length. From a fit of the lattice data to the operator product expansion, the 
quark condensate and the mixed quark-gluon condensate can be extracted. 
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1. Introduction 

Gauge invariant field correlators can serve as interesting examples for studying non- 
perturbative aspects of QCD in the framework of lattice QCD [1,2] or QCD sum 
rules [3-5]. In addition, they are a natural extension of the local condensates which 
appear in the sum rule approach, due to the use of the operator product expansion 
(OPE) [6]. Phenomenological implications of nonlocal condensates have previously 
been discussed in the literature [7-13]. The two most fundamental correlators are 
the gauge invariant field strength correlator and the gauge invariant quark correlator. 
The gluon field strength correlator has already been investigated in lattice QCD 
[14-16] and QCD sum rules [17,18]. These studies allowed to extract the correlation 
length of the gluon field. 

On the other hand, the gauge invariant quark correlator has so far only been 
measured on the lattice [19,20]. In this work, we therefore first present a QCD 
sum rule analysis of the gauge invariant quark correlator. To this end, the quark 
correlator is related to a heavy-hght meson correlator [21, 22] in the heavy quark 
effective theory (HQET) [23]. Thus the sum rule analysis will be related to previous 
investigations of the same HQET correlators [24, 25] in the hmit of the heavy quark 
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mass going to infinity. The correlation length of the quark field is then given as the 
inverse of the binding energy of the light quark. 

In the second, and more innovative, part of the present work, we perform a 
direct comparison of the lattice data [19,20] with the representations of the correla- 
tion function in the QCD sum rule approach. From the so-called phenomenological 
parametrisation — a single resonance plus the perturbative continuum as the simplest 
Ansatz — we again extract the correlation length of the quark field, in agreement 
with the direct sum rule determination. Fitting the lattice data with the theoretical 
correlator in the framework of the operator product expansion, we obtain estimates of 
the quark and the mixed quark-gluon condensate which are in reasonable agreement 
with sum rule phenomenology. 

Our paper is organised as follows: In the next section we discuss the relation 
of the gauge invariant quark correlator and the corresponding heavy quark current 
correlator. In section 3 we set up the expressions needed for the sum rule analysis 
and in section 4 we shall present our numerical results. Section 5 compares our 
results with recent lattice determinations of the quark correlator. In particular, a 
direct comparison of the OPE with the lattice data will be performed. Finally, in 
section 6, we conclude with a summary and an outlook. 

2. The gauge invariant quark correlator 

The central object of our investigation is the gauge invariant two-point correlation 
function of quark fields, 

V^(,{z) ^ {Q\T{<fMr\y.x)ql{x)m, (2.1) 

where the string operator (^"'^{y, x) ensures the gauge invariance of the correlator and 
is represented by 

with z — y — X. are the generators of SU(3) in the fundamental representation 
and V denotes path ordering of the exponential. In general, the gauge invariant field 
strength correlator could be defined using an arbitrary gauge string connecting the 
end points x and y, but in this work we have restricted ourselves to a straight line, for 
it is only in this case that the correlator (2.1) is related to a heavy meson correlator 
in HQET. 

Using the path-integral formalism, we are able to derive a relation between the 
correlator Vapiz) and the correlator of a local, gauge invariant current composed of 
a hght quark field q%{x) and an infinitely heavy quark field h1^{x) [22]. Analogously 
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to HQET, the heavy quark field h'^{x) is constructed from the field Q'^{x) with a 
finite mass mq in the limit mg — > oo, 

^«(^) = 4™oo ^ (1+5^)°/^ W ' (2-3) 

where is the four-velocity of the heavy quark. In the case of infinitely heavy 
quarks, the coupling to the gauge fields is given by the effective HQET action S'e// = 
/ dtkichufD^h with = 9^ — igAij_ [23]. In the free field case the heavy quark 
propagator is given by 

= (2.4) 

where is the zero-component of the velocity. In addition, we have the relation 
'^ti — with 1 2; I = \J — ie which follows immediately from the constraints 

z = z^v/v^ and v"^ — 1. In order to obtain a solution for the interacting case, we 
have to find the inverse of the operator iv'^D^. Using the following crucial relation 
obeyed by the string operator, 

yi'dl (Piy, x) = v^^ igT^A'l (y) 0(y, x) , (2.5) 

together with equation (2.4) one can show that the solution for the interacting case 
is found to be [22] 

{0\T{hl{y) h%x) e'^^ff]\Q) = S^p{z) (0|[Pe*^^''/o '^^-"^?(^+^^) . (2.6) 

The physical interpretation of this result is an infinitely heavy quark moving from 
point X to point y with a four-velocity acquiring a phase proportional to the 
path-ordered exponential. The limit mq — > 00 is necessary in order to constrain the 
heavy quark to a straight line and to decouple the spin interactions. 

The result (2.6) can be employed to relate the correlator (2.1) to correlators of 
gauge invariant local currents in HQET. To this end, we define the pseudoscalar and 
scalar heavy meson currents as 

jp{x) = ql{x){t^,Uh''p{x) and 3 s{x) = ql{x) K{x) , (2.7) 

as well as the correlators ^r(-z) with F = P or 5": 

Vr{z)^{Q\T{jr{y)Mx)m. (2.8) 

The two correlators could have also been defined with vector and axialvector currents, 
but as it was shown in ref. [24], the vector and axialvector correlators are equal to 
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the pscudoscalar and scalar correlators respectively in the heavy mass limit. This 
result is related to the fact that in the heavy mass limit the corresponding physical 
states become degenerate. Integrating out the heavy quark fields in (2.8), we obtain 
the expressions 

^p(^) = - \m\y){^-i^)r\y,x)q\xm s{z) , 

-^siz) = I {0\r{y)il+M''\y,^)Q\^m S{z) , (2.9) 

displaying the relation of the heavy meson correlators to the gauge invariant quark 
correlator of equation (2.1). 

In energy space the HQET correlator Vriw) obeys the spectral representation 

f 2 °° / \ \ 

^r(^) - E + / d>^ ' (2-10) 

k 



{Er,k -w-ie) J {X - w - ie) ' 



where w = vq and the Fourier transform of the correlators 'D-piz) in coordinate space 
is given by 

Vr{w) = ij d^z e"i' Vr{z) . (2.11) 

Er,k represents the energy of the bound states and fr,k is the coupling of the state k 
with quantum numbers F to the vacuum, 

{0\jrmHr,k)^fr,k. (2.12) 

The spectral densities are defined by pr{X) = l/irlmVriX + ie), and finally Wq 
is the threshold energy of the continuum contribution. Transforming the spectral 
representation (2.10) back to coordinate space, one obtains 

^^^^^ = (^""^"'^'^^'"^ ^ S{z)Vr{z) 

oo 

= S{z)[Y.fr,k^~'^''''^'^ + / ^^Apr(A)e-^l^l } . (2.13) 

K 

The factorisation of the heavy quark propagator and the relations (2.9) ensure a 
representation of the gauge invariant quark correlators T^riz) as given by the expres- 
sion inside the curly brackets. Let us already remark that the correlator decays as a 
simple exponential in the Euclidean region. The correlation length will therefore be 
given by the inverse of the lowest lying bound state energy. 
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3. The sum rules 



We now turn to the theoretical side of the sum rules which arises from calculating the 
correlator of equation (2.8) in the framework of the operator product expansion [3,6]. 
The perturbative contributions for the spectral density up to next-to-leading order 
in the strong couphng constant have been calculated in [24]. Up to C(m^) in the 
light quark mass, they are found to be: 



Pr(w) 



w 



n=0 



w J 



N, 



87r2 



4 + Ci.a(l7+|7r2-6L) 



(3.1) 



± wm 



4 + CFa (24 + |7r2-12L) - 2 + CFa(3 - 9L) | 



where a = as/n and L = \n{2w/fi). Equation (3.1) also defines the coefficients 
(ijjj which will be utilised below. As above and in all the following, the upper sign 
corresponds to the pseudoscalar and the lower sign to the scalar current. 

In order to suppress contributions in the dispersion integral coming from higher 
exited states and from higher dimensional operators, it is convenient to apply a 
Borel transformation Bt with T being the Borel variable. The Borel transformation 
also removes the subtraction constants which are present in the dispersion relation 
satisfied by the correlators. Some useful formulae for the Borel transformation are 
collected in the appendix. For the phenomenological side of the sum rules, equation 
(2.10), we only take the lowest lying resonance and approximate the spectral density 
by the perturbative expression (3.1), assuming quark-hadron duality for w > Wq. 
We then obtain 



Vt{T) 



+ 



j dAp^*(A)e-^/^. 



(3.2) 



Let us remark that equation (3.2) takes exactly the same form as the expression of 
(2.13) inside the curly brackets, which represents the gauge invariant quark correla- 
tor, if we identify 1/T with the Eucfidean space-time coordinate. 

The perturbative contribution relevant for the sum rules is the full correlator 
minus the corresponding continuum contribution: 

oo oo 

Vl\T) - V'j?{T, w^) ^ JdX ffr'W e'^'^ - J dX (A) e"^/^ , (3.3) 



where T>^{T, w^) denotes the perturbative continuum part. After the Borel transfor- 
mation the correlators satisfy homogeneous renormalisation group equations. Thus 
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we can improve the perturbative expressions by resumming the logarithmic contri- 
butions. With the help of the following general integral formula [26], 



dX A"-i In" — e-^/^ = T« (t) In' — 



Qn-k / ^rN 



(3.4) 



the perturbative contribution to the sum rules is found to be: 



(3.5) 



^(/)(3-n,7/) 



n=0 



dJjQ + a ( dF 



nl 



0'(3-n,|/) ^ 
0(3-n,y)^ 



^m(2T)^ 



where 0(q;, |/) = r(Q;) — r(Q;, y), 4>'{a, y) = d/da (j){a, y) and y — Wq/T. Some explicit 
expressions for the incomplete F- function r{a,y) and the functions (f){a,y), (f)'{a,y) 
can also be found in the appendix. In our notation /5i = (ll-/Vc — '^^f)l^ is the first 
coefficient of the QCD /3-function. The anomalous dimension 71 of both currents 
is easily found by reexpanding the running coupling and mass in terms of a{^) and 
■m{ii). The resulting expression is 



7i = 



nlL 



(3.6) 



where the first coefficient of the mass anomalous dimension, 71m = SCp/^,, has been 
used. The //-dependence of the correlators is due to the fact that the pseudoscalar 
and scalar currents are not renormalisation group invariant quantities. However, 
the product of m times the current, m-jr, is renormalisation group invariant in the 
full theory. Taking into account the additional multiplicative renormalisation of the 
heavy quark current in HQET [24,27], one again finds 71 = jim- 

Essential for QCD sum rule analyses are the non-perturbativc contributions com- 
ing from vacuum condensates. In the operator product expansion, the correlation 
function is expanded in powers of 1/w corresponding to higher and higher dimen- 
sional condensates. In our case the dimension three condensate (hh) vanishes since 
the heavy quark mass is infinite. After Borel transformation and up to operators of 
dimension five, the non-perturbative contribution takes the form [24]: 



m 3 ^ 



± 



32T^ 



(3.7) 



The next condensate contribution would be of dimension six. We have checked ex- 
plicitly in our numerical analysis that this contribution to the sum rule is small. 
Thus we shall neglect all condensate contributions for dimensions higher than five 
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in this work. For consistency, wc have also omitted the known order 0{a'^) contri- 
bution to the quark condensate [24], because the corresponding corrections to the 
perturbative part have not yet been calculated. We have, however, verified that also 
this correction only has a minor impact on our numerical results. 

In the case of the non-perturbative part, the scale dependence of the correlator 
is implicit in the /^-dependence of the quark condensate. Indeed, again m{qq) is 
scale independent and therefore (qq) scales inversely hke the quark mass yielding the 
same /x-dependence as for the perturbative contribution. For the mixed quark-gluon 
condensate, the next-to-leading order corrections have not been calculated and thus 
in the numerical analysis below, we shall neglect its scale dependence. 



4. Numerical analysis 

After equating the phenomenological and the theoretical contributions to the corre- 
lation functions, we end up with the following sum rule: 

KriT) = /2 e-^^/^ = Vp{T) - V'j?{T, w^) + V'^^T) . (4.1) 

The binding energy Ey can be obtained by dividing the derivative of the sum rule 
(4.1) with respect to — l/T by the original sum rule: 

(Pf (T) - (T, wl) + V^'^iT)) 

Analogously to equation (3.5) an expression for the derivative of the perturbative 
contribution can be obtained with the help of formula (3.4): 

d 



d{l/T) 



V'^{T)-V^{T,wl)) = (4.3) 



71 //9i 



The corresponding derivative of the non-perturbative contribution is easily calcu- 
lated from equation (3.7). Because the renormalisation of the correlators is multi- 
plicative, it is clear that the binding energy Ey is a physical quantity in the sense 
that it is scale and scheme independent. On the contrary, this is not the case for the 
decay constant /r, as we shall discuss in more detail below. 

Let us first consider the pseudoscalar correlator for a vanishing light quark mass 
niq — 0. As the central values for our input parameters we use (gg)(l GeV) = — (235± 
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Figure 1: Pseudoscalar binding energy as a function of l/T for different values of the 
continuum threshold Wq . 

20MeV)''^ for the quark condensate, {gqaGq) = (qq) with = (0.8 d= 0.2) GeV^ 
[28-30] for the mixed condensate and = 340 MeV. The latter value corresponds 
to three light quark flavours and as{Mz) = 0.119. In principle, the coupling constant 
in the next-to-leading order term could be evaluated at any scale fi. As our central 
value in the numerical analysis we have chosen fi = 1.2 GeV but we shall investigate 
the variation of fi below. In figure 1, we display the pseudoscalar binding energy 
Ep as a function of 1/T for different values of the continuum threshold Wq. A sum 
rule window, were both the continuum as well as the condensate contributions are 
reasonably small, can be found around 1/T = 2.5 - 3.0GeV"\ As can be seen from 
figure 1, best stabiUty is achieved for Wq ^ 0.9 GeV. Estimating the value of Ep in 
the given range 0.7 GeV < Wq < 1.1 GeV, we obtain: 

Ep = 450 ± 100 MeV (m^ = 0) . (4.4) 

Although the next-to-leading order QGD corrections to the sum rule of equa- 
tion (4.1) are very large, of the order of 100%, in the ratio of equation (4.3) they 
cancel to a large extent. To investigate the sensitivity of our result to higher order 
corrections, we next study the dependence on the renormalisation scale /i. The de- 
pendence of Ep on the renormahsation scale is shown in figure 2 for — 0.9 GeV 
and n = 0.8 GeV, 1.2 GeV and 2.4 GeV. One should not take n smaller than 0.8 GeV 
because then also the radiative corrections to the quark condensate become unac- 
ceptably large. We observe that inspite of the huge corrections to the correlation 
function, the scale dependence of Ep is relatively mild. Adding this uncertainty to 
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Figure 2: Pseudoscalar binding energy as a function of 1/T for Wq = 0.9 GeV and different 
values of the renormalisation scale fi. 



the error on Ep, our central result for Ep is found to be 

Ep = 450 ± 150 MeV {nig = 0) . 



(4.5) 



The additional uncertainty coming from the errors on the input parameters is small 
compared to the estimated uncertainty and can be neglected. 
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Figure 3: Pseudoscalar binding energy as a function of 1/T for Wq = 0.9 GeV, = 
1.2 GeV and different values of the light quark mass ruq. The absolute value of {qq) has 
been reduced in accord with the quark mass. 
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Let us now investigate the influence of a flnite light quark mass m^. In flgure 3, 
we have plotted Ep as a function of 1/T for Wq = 0.9 GeV, n = 1.2 GeV and three 
values of the light quark mass rUg, namely rriq = OMeV, 75MeV and 150 Me V. The 
last value is in the region of the mass of the strange quark and the intermediate 
value is interesting for comparison with lattice QCD results. It is well known from 
QCD sum rules that at the mass of the strange quark, the absolute value of the 
corresponding quark condensate is reduced by roughly 30% [4]. This reduction has 
been applied for obtaining the dashed curve in figure 3. Lacking further information, 
for the dotted curve corresponding to niq — 75 MeV the quark condensate has been 
reduced by 10%. 

Qualitatively, we find that the binding energy increases with increasing light 
quark mass and decreasing quark condensate. For niq — 150 MeV this increase turns 
out to be of the order of 50 MeV. Such a value is only about half of the mass 
splitting of B and Bg mesons in the S-meson system as well as D and Dg mesons in 
the D-meson system which is of the order of 100 MeV. Further comparison with the 
heavy meson systems will be presented below. The instability for 1/T ^ 3.5 GeV~^ 
results from the running quark mass mq{2T) which explodes in this region due to 
uncontrollably large higher order corrections. 

1.40 
1.20 

> 

O 100 
LU 

0.80 



0.60 

1.0 1.5 2.0 2.5 3.0 3.5 4.0 

1/T [GeV"'] 

Figure 4: Scalar binding energy as a function of l/T for different values of the continuum 
threshold Wq. 

Let us next turn to the scalar correlator. In figure 4 the scalar binding energy 
is shown for three different values of the continuum threshold Wq. We observe that 
generally the stability of the scalar sum rule as not as good as in the pseudoscalar 
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case. Here, the region of best stabihty is around 1/T = 2.0 GeV ^. Adding the 
uncertainty from the scale dependence, the scalar binding energy is found to be: 

= 1.0±0.2GeV (mg^O). (4.6) 

Compared to the pseudoscalar correlator, the dependence of Es on the light quark 
mass rrig is somewhat stronger. Approximately, we obtain Es{mq) ~ Es{Q) +mq for 
niq < 150 MeV. However, in the scalar case for increasing quark mass and decreasing 
quark condensate, the sum rule becomes less stable and thus we refrain from making 
more quantitative statements. 

To conclude this section, let us comment on the decay constants /p. As has been 
already discussed above, because of the renormalisation of the heavy meson current, 
/r depends on the renormalisation scale and scheme. Therefore, it is convenient to 
define rcnormalisation-group invariant decay constants /r- At the next-to-leading 
order, /r takes the form [24,25]: 

/r = /ra.(/.)-/^^ni + i^a) with ^ = ^ - ^^^^ + ^ • (4-7) 

As central values, from the sum rule of equation (4.1) we then find fp — 0.3 GeV^^^ 
and fs = 0.5 GeV^/^. Since the next-to-leading QCD correction to the correlators 
is very large, this is also true for the uncertainty on /r. Thus we shall not dicuss the 
heavy meson decay constants further. 

5. Comparison with lattice results 

In ref. [19] the two following quark-antiquark nonlocal condensates have been deter- 
mined on the lattice: 

Cv{x) = - ^ E mQimi'E)a,m x)4{x)]) . (5.1) 
/=i 

The trace in (5.1) is taken with respect to the colour indices and / is a flavour 
index for the light quarks. 7^ are Euclidean Dirac matrices, defined as: 7^4 = 7*^, 
7Ej = —iY {i — 1, 2, 3). Other insertions of gamma matrices than those containing 
1 and 7^ vanish by P and T invariance. 

In order to avoid confusion with the correlators introduced in section 2 the cor- 
relators Co{x) and C„(x) are called respectively "spin-zero" nonlocal condensate and 
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"longitudinal-vector" nonlocal condensate. They are simply related to the correla- 
tors Vr of equation (2.9) if the latter are continued to Euclidean space time: 



Vpix) = ^ (C,(x) + Coix)) Six) , 

-^six) = ^ (C^ix) - Coix)) Six) , (5.2) 



where Nf is the number of quark flavours. 

The lattice computations of ref. [19] have been performed both in the quenched 
approximation and in full QCD using the SU(3) Wilson action for the pure-gauge 
sector and four degenerate flavours of staggered fermions, so that the sum over the 
flavour index / goes from 1 to 4. In full QCD the nonlocal condensates have been 
measured on a 16^ x 24 lattice ait P — 5.35 and two different values of the quark 
mass: a • ruq — 0.01 and a • ruq — 0.02 (a being the lattice spacing). For the quenched 
case the measurements have been performed on a 16^ lattice at /3 = 6.00, using 
valence quark masses a ■ nig — 0.01, 0.05, 0.10, and at (3 — 5.91 with a quark mass 
a-ruq — 0.02. Further details, as well as a remark about the rehabihty of the results 
obtained for the longitudinal-vector nonlocal condensate, can be found in [19]. The 
scalar functions Cq and Cy introduced in ref. [19] have a more complicated spectral 
representation than Vp and Vs, since they receive contributions both from scalar and 
pseudoscalar intermediate states. Nevertheless, the correlator Cq has the advantage 
of not receiving perturbative contributions in the chiral limit — > with nig being 
the mass of the light quark: 

hm Coix) n n and lim C^(a;) ^ — , (5.3) 

where Nc is the number of colours. 

Since the available lattice results for C„(x) are not conclusive, in what follows 
we shall concentrate on the spin-zero nonlocal condensate Co(x). In ref. [19] a best 
fit to the data with the following function has been performed: 

Coix) = Aoexpi-iiox) + — . (5.4) 

The perturbative-like term Bo/x"^ takes the form obtained in the leading order in 
perturbation theory, if the chiral limit nig (see equation (5.3)) is approached. 
Using the ansatz (5.4), the correlation length Aq = l//Xo of the spin-zero quark- 
antiquark nonlocal condensate can be extracted. At the lightest quark mass a ■ nig — 
0.01 in full QCD one obtains the value Aq = 0.63io;i3 fm [19] corresponding to 
jjLo — 310 ± 80 MeV. Within errors this value agrees with the value for Ep obtained 
from the sum rule analysis. 
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Some observations are, however, necessary at this point. One should point out 
that the parametrisation (5.4) is a sort of "hybrid" parametrisation, since it contains 
a "perturbative" term Bq/x^, which should reproduce the behaviour predicted by per- 
turbation theory at short distances, and a "non-perturbative" term exp(— /iga:), 
which should reproduce the predicted exponential behaviour at large distances. A 
simple exponential term is not dominant in the range of physical distances where 
lattice data are taken, i.e., a; 0.2 0.8 fm. This could shed some doubt on the 
identification of /xq = I/Aq with the binding energy Ep, as determined in section 4. 

Nevertheless, at the distances where lattice data are taken the operator product 
expansion should still be a reasonable approximation. One can therefore try to fit di- 
rectly the OPE expression to the lattice data. In this case however a running coupling 
would lead to a Landau pole around 1 fm and one should confine the comparison 
to distances small compared to this scale. On the other hand most treatments of 
non-perturbative QCD are based on the assumption that at large distances only 
non-perturbative effects prevail. This can be achieved by freezing the coupling at a 
certain value for distances larger than a critical one. Given the fact that we have 
only few data points we consider here only the leading order terms in the OPE. In 
addition, the correlator is scheme dependent and to perform a consistent comparison 
between the lattice data and the OPE at the next-to-leading order, a perturbative 
calculation in the lattice scheme would be necessary. 

At the leading order the spin-zero correlator Co{x) up to operators of dimension 
seven is given by: 

Co{x)/Nf = K,{m,x) - I 1 + I my \ {qq) + — {gqaGq) , (5.5) 



TT^X 



where in the perturbative part through the Bessel function Ki{z) we have kept all 
orders in the quark mass. Fits for the quark and the mixed condensate from the 
different sets of lattice data of ref. [19] are given in table 1. The presented errors 
just correspond to the statistical la errors resulting from the fit if the x^/d.o.f. is 
normalised to one. is the lattice mass converted to physical units and ruq the 
mass appearing in the OPE of equation (5.5). Since we work at the leading order, 
the scale and scheme dependence of the quark mass are not under control. 

We note that the condensates come out with the correct order of magnitude. 
Qualitatively also the decrease of the quark condensate with increasing quark mass is 
found albeit somewhat stronger than expected from phenomenology [4] . A direct ex- 
traction of the quark condensate from the uncooled values of the spin-zero quark cor- 
relator [19] at zero distance led to the value {qq){1 GeV) = — (235 ± 15 MeV)^ in per- 
fect agreement with phenomenological determinations. We have again parametrised 
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P 


a ■ m 


niL [MeV] 


rriq [MeV] 


— {qq}' [MeVJ 




6.00, q 


0.01 


19 


38 ± 4 


244 ± 4 


0.47 ± 0.03 


5.91, q 


0.02 


33 


73 ± 11 


198 ± 12 


0.44 ± 0.06 


o.OO, q 


A A r 

0.05 


A/^ 

96 


1 TA 1 1 A 

150 ± 14 


186 ± 15 


A O 1 AAA 

0.68 ± 0.09 


6.00, q 


0.10 


192 


276 ± 22 


175 ± 20 


0.87 ±0.12 


5.35, f 


0.01 


20 


28 ±2 


168 ±4 


0.47 ±0.04 


5.35, f 


0.02 


33 


54 ±6 


177 ±9 


0.45 ± 0.04 



Table 1: Determination of the condensates by direet eomparison of lattice data with the 
OPE. q denotes quenched approximation, f fuh QCD with 4 hght fermions, is the input 
lattice mass and nig the OPE mass appearing in equation (5.5). 



the mixed quark-gluon condensate through (gqcrGq) = MQ{qq). The dependence 
of Mq on the quark mass is controversial in the hterature [31,32]. Prom our fits 
we obtain an increase of Mq with increasing quark mass, supporting the findings of 
ref. [31]. 

The non-perturbative part of the OPE (5.5) can be expressed at short distances 
by the Gaussian — (^g') exp(— MqX^/16) which at large distances displays an expo- 
nential falloff. Pitting the lattice data for the quenched calculation with a-ruq — 0.01 
to the perturbative part of (5.5) and the Gaussian non-perturbative contribution, 
we find ruq = 33 ± 2MeV, -{qqY'^ = 254 ± 2MeV and = 0.74 ± 0.02 GeV^. 
This value for Mq is more compatible with phenomenological determinations [28-30] . 
However, the result shows that higher order corrections in the OPE have some im- 
portance and it gives an indication about the systematic uncertainties. 

Another way to improve the hybrid expansion (5.4) consists in taking the higher 
resonances into account in the same way as it is done in the sum rule technique 
as developed by SVZ [3], namely by including the perturbative continuum above a 
threshold Wq. In leading order for the spin-zero nonlocal condensate this yields: 

((2 + 2w^x + {w^xf)e-<^ - (2 + 2w^x + {w^xf)e-^'^'') 

+ ((1 ± w^x)e-<^ + (1 + w^x)e-^o-^ . (5.6) 

In figure 5 we have displayed the lattice data for the quenched calculation with 
a • niq — 0.01. The dashed curve is the fit with the OPE (5.5) and the parameters 
given in the first row of table 1. The solid curve is a fit with equation (5.6) where for 
simplicity, and lacking enough data points, we have neglected the scalar resonance 
by setting fs — and Wq — 0. This means that in the scalar channel only the 
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0.12 I • • • • < • • • • < • • • • < • • • r 




0.02 I — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — ' — 

1.0 2.0 3.0 4.0 5.0 

X [GeV"'] 

Figure 5: Lattice data from [19] with a ■ ruq = 0.01, quenched calculation. Dashed 
curve: OPE with condensates of table 1, first row, solid curve: equation (5.6) with fp = 
0.245 GeV^/^ Ep = 411 MeV, = 1.0 GeV, fs = and = 0. 

perturbative term was taken into account and no resonance was singled out. Using 
in addition rriq = 38 MeV, the fit leads to fp = 0.245 ± 0.002 GeV^/^ 

£;p = 411±3MeV, (5.7) 

and Wq = 1.0 ± 0.01 GeV. Within the uncertainties, these results are in good agree- 
ment to the sum rule determination of the last section. 

In other lattice investigations [33,34] of the heavy meson systems the main inter- 
est was to determine the leptonic decay constant of the heavy light meson in HQET. 
In both cases it was found that in order to isolate the ground state contribution one 
has in some way to suppress the higher state contributions. This was done by two 
different smearing mechanisms. Though the authors claimed that on the lattice the 
mass gap Ep is not a physical quantity and therefore divergent in the continuum 
limit they determined this quantity for their lattice spacing. The difference of the 
mass gaps Es — Ep however is a physical quantity and must therefore have a con- 
tinuum limit. In table 2 we collect the results for these quantities for the different 
approaches and compare to our findings presented above. 

6. Conclusions 

In this paper we have investigated the gauge invariant quark correlator. We have set 
up a relation between the correlator and a corresponding correlator of gauge invariant 
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Ker. 


P 


a [GeVJ 




hs — -bp [GcVJ 


To A 1 

[34] 


5.74 


1.53 


1.93^ 


0.69 


[34] 


6.0 


2.03 


1.58* 


0.49 


Fool 

[33] 


6.U 


2 


1.3 




[19] 


6.0 


1.9 


0.31 




this work L 


6.0 


1.9 


0.41 




this work SR 






0.45 


0.55 



Table 2: Pscudoscalar mass Ep and difference between scalar and pseudoscalar mass in 
different approaclies for different values of the strong coupling /3 = 2Nc/g1 and lattice 
spacing a. This work L is an analysis of the data lattice data from [19] according to 
equation (5.6), this work SR the sum rule analysis of section 4. In the lattice calculations 
of [33, 34] the pseudoscalar mass is not defined for a 0. 

currents in the heavy quark effective theory. The relevant currents are interpolating 
fields of pseudoscalar and scalar heavy-light mesons. With the method of QCD sum 
rules and cooled lattice data it is possible to extract the mass gap between the heavy 
quark pole mass and the lightest bound states of each current. This mass gap is a 
physical quantity in the sense that it is scale and scheme independent in perturbation 
theory to all orders in the strong coupling constant. The resulting values have been 
found to be Ep = 450 ± 150 MeV for the pseudoscalar state and Eg — 1.0 ± 0.2 GeV 
for the scalar state. The corresponding correlation lengths are: ap = 0.44to;ii fm and 
tts — 0.20lo;o3f™- For the cooled lattice data of [19,20] an analysis of a perturbative 
term plus an exponential gave the correlation length ap = 0.63lo;i3fm; the modified 
analysis with a single resonance plus continuum of equation (5.6) yields the value 
Ep = 411 MeV, corresponding to ap = 0.48 fm. 

These results can be compared to the spectrum of heavy pseudoscalar mesons. 
Let us make the very simple assumption that up to corrections of 1/ rriQ where mq is 
the heavy quark pole mass, the heavy meson mass is equal to the heavy quark mass 
plus binding energy: 

Ep- — = Mp-mq. (6.1) 
rriQ 

Here Mp is the mass of the pseudoscalar meson and c is a constant. Assuming in 
addition that this relation is vahd for both the B and D mesons, we can solve for 
Ep and c. Using nif, — 5.0 ± 0.2 GeV, rric — 1.8 ± 0.2 GeV and experimental values 
for Mb and Md, as central values for Ep and c we obtain: 

Ep = 400 MeV and c = 0.6 GeV^ (6.2) 

The result for Ep is in good agreement with the other determinations presented 
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above. Nevertheless, we should remark that the latter values are very sensitive 
to the heavy quark pole mass and with the estimated errors on mf, and rric, the 
uncertainty on Ep is of the order of 100%. In addition, from the value of c we see 
that the assumption of small l/mq corrections is not valid in the charm case. Still, we 
find the very qualitative agreement of our results with the spectrum of pseudoscalar 
mesons noteworthy. 

For lattice simulations of the gauge invariant quark correlator where the higher 
states were suppressed by a smearing procedure the mass gaps are said to diverge 
in the continuum limit. Though there is no indication of such a divergence with 
increasing /3 (see table 2) the very large values for Ep found in [33] and [34] are 
definitely not physical. The difference Es — Ep is however a physical quantity. It 
has been estimated in [34] and the values listed in table 2 agree with the result from 
the sum rules better than expected in view of the large errors. 

The lattice data for Co{x) have directly been compared to the operator prod- 
uct expansion in leading order QCD. The results for the continuum values of the 
quark mass, the quark and the mixed condensates are well compatible with the val- 
ues determined from other sources. The results from the quenched approximation 
are nearer to the continuum values than those from full QCD. The known decrease 
of the modulus of quark condensate with the quark mass is also confirmed by the 
lattice calculations. The ratio of the mixed to the quark condensate Mq comes out 
to be the same in the quenched and full QCD. There is some controversy on the 
dependence of Mq on the quark mass [31,32]. The lattice data support an increase 
of Mq with the quark mass which was found in [31]. For distances smaller than 
approximately 1 fm where the operator product expansion can still be applied a 
good fit to the non-perturbative part of the correlator is also given by the Caus- 
sian —{qq) exp(— MqX^/16) which has previously been used in the nonlocal sum rule 
approach [11-13]. 

In our opinion a direct comparison of lattice QCD simulations with QCD sum 
rules in the framework of the operator product expansion opens a novel route to 
augment our knowledge on low-lying hadronic states and non-perturbative QCD. 
Taking the results of this work as encouraging, we intend to further pursue this 
approach in the future. 
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A. Appendix 



The Borel transformation in HQET is defined as 

yn+l / n 



Bt = lim 



—w 



n 



> fixed . 



(A.l) 



o,n^oo r(n + 1) \dw J ' 
Using this definition, a central formula for the Borel transformation is found to be 

1 1 



Bn 



{E-w- ieY V{a)T<^-^ 



-E/T 



(A.2) 



Below, we have collected some analytic formulae for the incomplete Gamma 
function and the functions (f){a,y) and (f)'{a,y) defined as 



0(Q;,y) = r{a) - r{a,y) , 
d 

which are helpful for the numerical analysis of the sum rules. 



(A.3) 



my 

r(3,y 

r'{2,y 
r'(3,y 

(l)'{a,y 
0'(3,y 



= e-y{l + y) 

= T{0,y) + e-y[l + {l + y)\ny] 

= 2r(0, y) + e-y[3 + y + (2 + 2y + y^) \ny] 



= Tin) 1 



n-l yk ' 



E 

fc=0 



n = l,2,... 



—t j.a—1 



dt \nt e'H' 

= --fE - r(0,i/) - \ny 

= 1 - 7£ - r(0, y) - (! + (! + y) Iny) 

= 3 - - 2r(0, y) - e"^ (3 + y + (2 + 2y + y^) In y) 

= ll-67B-6r(0,y) 

-e-y (ll + 5y + y^ + (6 + 6y + Sy^ + y^) Iny) . 



(A.4) 
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